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A ROBUST PSEUDO-INVERSE SPECTRAL FILTER APPLIED TO THE EARTH 
RADIATION BUDGET EXPERIMENT (ERBE) SCANNING CHANNELS 


ABSTRACT 

Computer simulations of a least squares estimator operating on the ERBE 
scanning channels are discussed. The estimator is designed to minimize the 
errors produced by non-ideal spectral response to spectrally varying and 
uncertain radiant input. The three ERBE scanning channels cover, respec- 
tively, a shortwave band from 0.18 to 5 microns, a longwave band from 5 to 
50 microns, and a "total" band from 0.15 to beyond 1000 microns, from which 
the pseudo-inverse spectral filter estimates the radiance components in the 
shortwave band and a longwave band (5 to 1000 microns). 

The radiance estimator draws on instantaneous field of view (ITOV) scene- 
type information supplied by another algorithm of the ERBE software, and on a 
probabilistic models of the responses of the scanning channels to the 
IP0V scene types for given Sun-scene spacecraft geometry. 

The pseudo-inverse spectral filter has been found to be stable, tolerant 
of errors in scene identification and in channel response modeling, and, in 
the absence of such errors, to yield minimum variance and essentially unbiased 
radiance estimates. 
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INTRODUCTION 


The "radiation budget" is the solar irradiance absorbed in a given region 
of the Earth-atmosphere system minus the radiant emittance to space averaged 
over a specified time. Knowledge of how the radiation budget varies tempo- 
rally and spatially over the globe is fundamental to understanding the energy 
sources and sinks within the Earth-atmosphere system. 

Satellite observations of the radiation budget date from 1959, with the 
Explorer 7 providing the first. There followed TIROS 2, 3, 4, and 7; Nimbus 2- 
ESSA 3, 5, 7, and 9; Nimbus 3; ITOS 1; NOAA 1; and Nimbus 6 a^d 7. ihe Srth 
Radiation Budget Experiment (ERBE), benefiting from the experience gained by 
these earlier efforts, is expected to yield greatly improved measurement accu- 
racy through emphasis on instrument characterization and calibration, both 
pre-flight and in-flight. Furthermore, the three satellites of the ERBE 
system should reduce errors caused by aliasing of Earth flux field variations 
with sampling periodicities. 

The National Aeronautics and Space Administration (NASA) is planning the 
ERBE mission for the mid-1 980's to provide an initial research data set and to 
demonstrate the feasibility of routine monitoring of radiation budget param- 
eters for use in climate analyses and projections. The ERBE radiometers will 
be flown on two of NOAA’s Sun-synchronous TIROS N satellites and a NASA- 

developed mid-inclined Earth Radiation Budget Satellite (ERBS) to be launched 
from the Space Shuttle. 

The three satellites are to have nearly identical ERBE instrumentation: 

( 1 ) A total wave band solar monitor 

(2) Limb-to-limb Earth viewing radiometers, one for the shortwave band 

(0.18 to 5 microns) and one for the total band 

(3) Medium field of view (10° Earth central angle) nadir viewing 

radiometers with the same bands as the limb-to-limb sensors 

(4) Three Earth scanning radiometers covering, respectively , a shortwave 

band (0.18 to 5 microns), a longwave band (5 to 50 microns), and a 

"total" band (0.15 to 1000+ microns)* 

Ihe principal mission of the scanners is to monitor emitted and solar 
reflected radiant exitance from 250- to 1000-kilometer regions. Any use of 
these data sets to improve understanding of the energy sources and sinks and 
their action on the atmosphere and oceans could conceivably be misled not only 
by the magnitude of the errors in the data, but also, a fortiori, any correla- 
tions the errors may have with the processes under study. The scanning chan- 
nels do not respond uniformly to all wavelengths within their respective wave- 
bands; thus, to the extent that variable spectral response is not accounted 
for, there will be scanner measurement errors correlated with spectral changes 
in the scenes in view— notably, the spectral changes associated with changes 
in cloudiness, snow cover, vegetation, and solar zenith angle. 
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Some of the planned ERBE scientific investigations which could be susceptible 
to variable scanner spectral response are: 

(1) Albedo-solar zenith angle dependence as a seasonal feedback mechanism 

in climate models 

(2) Effect of cloudiness on radiation budget 

(3) Radiation budget-atmospheric circulation correlations 

(4) Radiation parameterization in global circulation models (GCM's) 

(5) Radiation budget as diagnostic for GCM's. 

In this paper, we present a technique aimed at (1) reducing scanner 
spectral response effects and (2) quantifying the remaining errors of this 
kind in terms of ejq>licit models of the scanner channels' spectral response 
and stochastic models of the Earth spectral fields. We feel that making the 
models and assumptions explicit is necessary for the intelligent use of the 
data products, and also allows updating the technique whenever it is necessary 
or desirable to alter the models because of changes in the instrument charac- 
teristics or improvement in our knowledge of the Earth spectra. 


BACKGROUND 

Description of the ERBE Scanner 

The scanning channels are configured as three parallel telescopes which 
scan the Earth from limb to limb. The lines of sight are confined to a plane 
perpendicular to the spacecraft (nominal circular orbit) velocity vector. Ebr 
special, occasional studies, the scan plane is rotated about the local verti- 
cal to other orientations for brief periods. In the principal or nominal scan 
mode, the scan period is 4 seconds, which assures sampling overlap from scan 
line to scan line. Sampling overlap along the scan line is also realized by a 
sampling rate per channel of 30 hertz. The instantaneous fields of view for 
all scanning channels measure, at the spacecraft, 3° along scan and 4.5° along 
ground track. The fields of view (clipped diamonds) project to the Earth at 
nadir as nearly coincident areas of dimensions 47 by 37 km for the ERBS mid- 
inclination satellite, and 65 by 51 km for the TIROS satellites. 

The telescopes employ Cassegrain optics with polished aluminum primary 
and secondary mirrors. The shortwave channel has a filter and a window of 
Suprasil-W fused quartz, one-half millimeter and one millimeter thick, respec- 
tively. The longwave channel has a diamond filter. The total channel, 
although not filtered, still suffers appreciable spectral nonuniformity by 
double reflection off the mirrors. The current preliminary models of the 
spectral responses of the three channels are shown in figures l.a-l.f. The 
spectral absorption of the thermistor flake detector is factored into the 
response curves. 

Fortunately, inferrence of the longwave radiance is not solely dependent 
on the longwave channel; the total and shortwave channels in combination yield 
good longwave radiance estimates, in part because these two channels have very 


3 



similar spectra curve shapes in the near ultraviolet and the visible bands. 
This spectral similarity derives from the nearly ideal spectral transmittance 
of Suprasil-W in the near ultraviolet and visible, and from the presence of 
aluminum mirrors in the optical trains of both the total and shortwave 
channels. 

The longwave channel filter was chosen to be spectrally similar to that 
of Nimbus 6 and 7 ERB longwave channels, so that more direct connection might 
be made between ERBE and ERB measurements than could be made with spectrally 
dissimilar channels. The longwave channel, as well as offering a measure of 
redundancy, is intended to assist in cloud detection and classification. The 
redundancy feature is significant; the shortwave channel cannot be calibrated 
flight against blackbody simulators as the longwave and total channels are, 
and Suprasil— W transmittance in the near ultraviolet is known to degrade under 
ultraviolet and energetic electron flux. (Preflight flux aging is expected to 
desensitize the Suprasil to orbital flux.) In-flight checks of the shortwave 
channel are made by observations of space, an incandescent lamp, and of the 
Sun’s reflection from the "Mirror Attenuator Mosaic" (MAM), which functions as 
® diffuser of solar flux. The MAM and the incandescent lamp are not consid- 
ered of the same order of accuracy and precision as the blackbody simulators, 
however. The total channel is also checked against space and the MAM, as well 
as a blackbody simulator. The longwave channel is in-flight calibrated 
against space and, as mentioned before, a blackbody simulator. 


The Scene Identification Algorithm 

Evidently, the uncertainties in estimating radiance inputs to channels 
with nonuniform spectral response cannot be bounded, much less minimized, 
unless the uncertainties in the spectral content of the inputs are bounded in 
some sense. The ERBE scene identification algorithm, developed for the pur- 
pose of correcting for the anisotropy of the Earth's outgoing radiation, is 
also applicable to the scanner spectral correction problem. 

The scene viewed by the scanner is classified, pixel by pixel, as water; 
land, but not desert; snow; desert; land-water mix; land-snow mix; water under 
partly cloudy sky; non-desert land under partly cloudy sky; desert under 
partly cloudy sky; land-water mix under partly cloudy sky; land-snow mix under 
partly cloudy sky; water under mostly cloudy sky; non-desert land under mostly 
cloudy sky; desert under mostly cloudy sky; land-water mix under mostly cloudy 
sky; land-snow mix under mostly cloudy sky; or overcast sky. The clear sky 
scene types, the first six in the above list, are considered to have cloud 
amounts from 0 to 5 percent, the partly cloudy scene types from 5 to 
50 percent, the mostly cloudy scene types from 50 to 95 percent, and the 
overcast condition more than 95 percent cloud amount. The mixed undersurfaces 
range from 33 to 67 percent to 67 to 33 percent composites of the two pure 
undersurface types. In general, undersurface types are defined to resolution 
elements 2.5° by 2.5° in latitude and longitude. Cloud amount, however, is 
determined pixel by pixel. 

The undersurface type distribution over the globe is defined a priori, 
while the cloud amount in the instantaneous field of view (IPOV) is con- 
structed from the scanner measurements. 
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As discussed in a following section, statistical populations of radiation 
spectra, generated by radiation transfer models, are associated with the scene 
types (sub-divided into polar, mid-latitude, and tropical regions and 
according to Sun-IFOV-spacecraft angles)* The pseudo-inverse spectral filter 
operates on numerically integrated channel responses to each spectrum and on 
statistical parameters of the spectral population to form optimal radiance 
estimates as a linear combination of the channel measurements. 

In view of the dependence of the spectral correction on the constructed 
IFOV scene type, it seems appropriate to describe the ERBE scene identifica- 
tion algorithm in some detail. 

The scene identification algorithm relies upon models which predict the 
shortwave and longwave radiances to be observed by ERBE from each scene type. 
The predictive models are constructed from archived Nimbus 7 ERB scanner data 
and THIR 11 -micron window channel observations. The window channel serves to 
estimate cloud cover fraction by means of an objective threshold technique, 
while the ERB scanner measurements accumulate statistics on the radiances 
associated with each scene type as functions of the exiting ray direction and, 
for shortwave, the incident solar ray direction (ref. 1). From the predictive 
model thus constructed, the ERBE scene identification algorithm computes the 
probability that each ERBE scanner IFOV overlays a "clear, 11 "partly cloudy," 
"mostly cloudy," and an "overcast" area. The undersurface is specified 
according to geography. 

Taylor and Stowe (ref. 1) have concluded that except for land areas at 
night (often classed as cloudy), the cloud amounts derived from THIR are 
accurate to 10 to 20 percent, and that the amount of low clouds is the most 
difficult to estimate accurately. 

As things stand at present, we must expect appreciable confusion between 
clear targets and those underneath warm (i.e., low contrast in the infrared) 
clouds. Plans are being made to use GOES visible and infrared measurements of 
ERB/THIR targets to improve the clear surface-warm cloud discrimination. The 
poor discrimination between clear land and clouds at night is not considered 
to seriously degrade either radiation anisotropy modeling or spectral response 
correction. 

Although the Nimbus 7 ERB shortwave and longwave scanning channels pass 
the same nominal wavebands as the ERBE shortwave and longwave scanners, there 
are some possibly significant differences in detail which, of course, tend to 
soften the credibility of ERBE scene identification built on ERB data. The 
ERB shortwave channel optical train includes two mirrors, a 2-millimeter-thick 
Suprasil-W filter, and a pyroelectric detector. The ERB longwave channel has 
two mirrors plus a third chopper/beam diverter mirror, a diamond filter very 
similar to the ERBE longwave filter, and a pyroelectric detector. The spectral 
responses of the Nimbus 7 scanning channels have not been measured. 
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THEORY OF THE SPECTRAL CORRECTION 


The spectral radiance L(X) incident on the ith scanning channel 
(i=1*2,3) produces a measurement given by 


Y. 

i 



r. (X) L(\) dX + e. 
x x 


( 1 ) 


where r^X) is the spectral response, X is the wavelength (|jm), and e- 
is the measurement error arising from all sources except uncertainty in the 
normalized spectral radiance. 



L ( X) dX 


(2) 


The measurement Y i characterizes the integral of spectral radiance 
L(X) over a waveband labeled "i, M defined by: 


i waveband ( \m ) 


channel/waveband name 


1 

2 

3 


0-5 
5-a> 
0 — 00 


shortwave 

longwave 

total 


The spectral response r^(\j is dimensionless and, for an ideal instrument, 
is unity within the channel waveband and zero elsewhere. Losses in the 
optical train of the actual instrument reduce the spectral response in the 
passband to a nonuniform value below unity (fig. 1), and shortwave filter and 
window transmission at 50 pm and beyond leaks longwave^ into the shortwave 
measurement (fig. 2). 

The shortwave channel is predominantly sensitive to reflected solar 
energy, and the longwave to Earth/atmosphere emitted energy. 

The shortwave radiance, 

f 5 

L = I L(\) d\ ( 3 ) 

Jo 

and the longwave radiance. 


^Cess, 


Robert D.: 


Private communication. 
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L(X) dX 


(4) 


are estimated by 
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i=c 


1,2; c 


1 , day IFOV 

2, night IFOV 


(5) 


where L 1 and L 2 are the estimates of L 1 and L 2 , and the matrix of B's 
is a variable to be determined for each IFOV. The estimate error is 



( 6 ) 


2 2 

We choose to minimize the expected value of + n 2 » or 

t'l-t [(Z^ M - -J 2 

j=1 j-1 


a necessary condition for which is 


3B 


nm 


= 0? n = 1,2f 


in ■" c , • • « , 3 


or, since the operators E (expectation) and 


3B 


commute. 


nm 


8B 


nm . . 
J=1 


2 

L 


= 0 


or, by use of equation (7), 


( 8 ) 


E 


[B 1c V J - L ,j l m • S 'n * [B 2i Y J ' Y m * S n 


= 0 


7 



where the 5's are Kronecker deltas; thus 



which are two sets of 4-c equations in 2x(4-c) unknown B’s expressing a 
necessary condition for minimum sum of squared error, v\ + r| 2 # For a given 
IFOV, the B's are constant, thus eq. (9) becomes 1 2 


E 



do) 


or, in matrix notation. 



e{y y t } 

li 

l« C 

e{l y! 
1 n 1 

where 

n = 1 or 2 

and 

Y and 5 are 

of the 

Y 1 s and 

B n -s. 

Therefore, 


( 11 ) 

the column vectors (dimensioned 4-c) 


5 n = [e{y y t }] _1 E {L n y} 


for nonsingular e{y Y T } . Let 



day IFOV 


night IFOV 


( 12 ) 


(13) 


8 


L 


(14) 




"" A. “ 

L 1 

A 

L 1 


? L = 

A 

L 0 


L 

L 2 - 


L 2 J 


Then 


[B] = [e{y Y T }] _1 e{y L T } 


and, from equation (5), 


A_ _ 

l t = y t (b] 


for nonsingular e{y Y T } . 

Equation (16) is the formal machinery for optimal radiance estimation 
the presence of all instrument errors. 

The integrated response A, of the ith channel to the normalized 
spectral radiance 1 ( X ) (eq. (2)) is 


A, = J r. (X) 1(A) dA; i = c,...,3 
•'0 


and so, 


Vi = V L , + L 2> + 


The elements of the matrix e{y L T } become 


e{y l 1 =b{[a (l, + lJ + e 1 l | 

1 m n J lL m 1 - 1 2 ' m J n ; 

= E ( A m L L n } + E { £ m E n }; m = c,...,3; n = 1,2 


where L - L 1 + Lj, the total radiance. The elements of the matrix e{y Y 
become 


(15) 


(16) 


in 


(17) 


(18) 


19) 

} 
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( 20 ) 


e{y Y } = e{(a L + e )(a L + e )} 
s m J lv s s m nr J 

2 

= e K A m L } + E {e A l} + e{e a l} + e {e e } 
sm 1 m s J 1 s m J l sm J 

S/m = c, . • . / 3 


We note that modeled instrument noise gets into the influence matrix [B] 
through the last term in equation (20), 

In outline, our approach to radiance estimation is to prepare the influ- 
ence matrices [B] (eq. (15)) in advance of the flight for use in equation (16)* 
In this way, demand on computer resources during operational data processing 
is minimized, and assurance of algorithm stability for all conditions is 
obtained in advance. 

In general, a number of theoretically constructed Earth spectral radi- 
ances and probabilistic weights are associated with each of the scene types 
(subdivided into polar, mid— latitude, and tropical regions, and according to 
Sun-IFOV-spacecraf t angles) in the repertoire of the ERBE scene identification 
^^- 9 °rithm. These fields and weights, with the scanning channels* spectral 
responses r.(\), determine the influence matrix [B] by equations (19), (20), 
and (15). 1 

Numerical instability of the algorithm is possible (but easily avoided) 
as a result of amplification of truncation errors by an ill-conditioned 

e{y Y } matrix in equation (15). The ill-condition occurs as the matrix 

e{y Y } approaches singularity. The following demonstrates that singularity 
occurs if and only if the measurement from any active channel is modeled as 
deterministically related to the other active channels* measurements. 

Let 


p = the number of model spectra in the IFOV 
P k = the model probability of the kth spectrum; k = 1,...,p 
— sth channel model measurement of kth spectrum 


x 1 
cn 

1 


Y S1 

X s2 

• 

ii 

Y s2 

. 

. 

. 

X 

sp 


1 

* w 

• . 
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Then, the elements of e{y Y T } are 


P 

{Y Y } = V [ p , Y Y 1 = X • X 
1 s m j 2^ L k sk mk J s m 

k=1 


The vectors are linearly dependent if and only if (ref. 2) the Gram's 

determinant 


X • X =0 
s m 


which, for three active channels, is equivalent to 


x 3k = (X *2k + ^ X 1k ; a 'P constant; k = 


or 


*3* ' “ Y 2k + f»,k 


( 21 ) 


and for two active channels (nighttime IPOV), 


Y 2k = “ Y 3k 


( 22 ) 


An ideal instrument, with flat spectral response and e i = 0, would always 
have a singular Gram's determinant. In general, a vanishing Gram's determi- 
nant implies channel redundancy. The ERBE scanner is not sufficiently close 
to the ideal for channel near redundancy to pose a serious problem, but 
measures are taken to guard against it, nevertheless. 


SPECTRAL RADIANCE MODELS 

Both the longwave and shortwave spectral radiance modeling use the 106 
radiosonde-measured atmosphere models taken from Wark et al. (ref. 3). The 
atmospheric models are profiles of pressure, temperature, water vapor, and 
ozone for all seasons and region types of the Northern Hemisphere, repre- 
senting a very wide range of meteorological conditions. 


Longwave Spectral Radiance Models 

The spectral band 5 to 200 |im (2000 to 50 cm is divided into intervals 
10 cm wide. Atmospheric transmittance for each interval was calculated 
using the quasi-random band model of absorption (ref. 4). It was shown by 
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Tiwan and Gupta (ref. 5) that spectral transmittances of several infrared 
bands of atmospheric gases computed using this band model compared extremely 
well with laboratory measured values. Spectral parameters of individual 
absorption lines for the computation were obtained from the AFCRL line param- 
eters compilation (ref. 6). Continuum absorption by water vapor and its tem- 
perature dependence have been taken into account as outlined by Roberts et al. 
(ref. 7). For spectral regions where bands of different constituents overlap, 
overall transmittances were obtained by multiplying the individual values. 
Absorption due to all major and minor constituents has been taken into 
account. Spectral radiance for overcast atmospheres is obtained by treating 
the cloud top as the underlying surface and considering only that part of the 
atmosphere which lies above the cloud. For partly cloudy conditions, the 
spectral radiance is obtained as the cloud-fraction weighted sum of the 
overcast and clear spectral radiances. 


The model atmosphere is divided into homogeneous layers with boundaries 
at 900, 800, 700, 600, 520, 450, 400, 350, 300, 250, 150, 100, 50, 20, and 
10 millibars, truncated as appropriate by the underlying surface/cloud height. 


The exiting spectral radiance is output at wavelengths 5 pm, from 6 
through 50 pm at 2-pm intervals, from 50 through 100 pm at 10-pm intervals, 
and from 100 through 200 pm at 20— pm intervals. 


The spectral radiance is output at zenith angles 0°, 30°, 45°, 60°, and 
75°. Spectral radiant exitance is also computed. 


Shortwave Spectral Radiance Models 

The upwelling spectral radiant exitance and spectral radiance at 50-km 
altitude, the assumed "top of the atmosphere," is computed by a method based 
on the finite difference method of Barkstrom (ref. 8) but generalized to 
include azimuthal variations, spectral variations for the entire solar spec- 
trum, and specular reflection at the surface (ref. 9). The atmosphere is 
approximated locally as a plane-parallel layer with detailed vertical profiles 
of the physical and optical properties of the important atmospheric gases and 
aerosols. The cloudy atmosphere case is approximated by a plane-parallel, 
homogeneous layer. 


The formulation of the problem is patterned after the "planetary atmos- 
phere problem" given by Chandrasekhar (ref. 10). The system is assumed to be 
(1) composed of a series of plane-parallel layers so that physical properties 
vary in the vertical direction only, (2) illuminated at the top by a colli- 
mated beam of unpolarized solar radiation, and (3) bounded by a lower surface 
with specified reflective properties. Two additional simplifying assumptions 
are made: Thermal emission can be neglected and, for a clear atmosphere with 

an aerosol haze layer, polarization effects can be neglected (ref. 11). 

A prohibitive computer cost would be incurred in generating spectral radi- 
ances from 106 model atmospheres with several different undersurface types for 
several solar zenith angles at several exit ray zenith angles and a few differ- 
ent exit ray azimuths. Some judicious choices must be made concerning which 
spectral radiant exitances are to be expanded into spectral radiances, and how 
to represent the remaining spectral radiances in terms of the expanded ones. 
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We are presently implementing and testing an objective procedure for making 
this decision. The procedure is founded on the (thus far tested) idea that 
"similarity" (in some sense) of two spectral radiant exitances bespeaks simi- 
larity of their expansion into spectral radiances. Some number of radiant 
exitances having minimum clustering in "similarity" would then be expanded 
into spectral radiances, which would best span the spectral radiance 
variations. 

Two candidate definitions of "similarity" are being tested. The first is 


(D f 5 pm / 

SI. =1 F. (X) • F . ( X) dX/f 

^ Jn 1 3 


where 


i 


( C S ym 2 / 

f 5 Vim ) 

= max jJ Fj(X) dX, J 

F (X) dX 

0 3 ) 


and F^, Fj are spectral radiant exitances; and the second is 

(2) r 5 pm / /r 5 m 2 

J 0 1 / Vo 


'(X) dX 


L 


5 ym 


FT ( X) dX 
D 


The first definition of "similarity" exacts a penalty for differences in the 
integrated radiant exitances, while the second does not. Both similarities 
range from zero to one. 

Once a suitable definition of "similarity" is discovered, some number of 
minimally clustered spectral radiant exitances for each solar zenith angle and 
scene type category will be expanded into spectral radiance angular distribu- 
tions. The remaining spectral radiances are to be represented by the short- 
wave and total channel responses as approximated by a linear interpolation 
between nearest neighbors in the "similarity" continuum. 

The notion of "minimum clustering" requires definition. The first two 
spectral radiant exitances to be expanded into spectral radiances are two with 
the least pair-wise similarity. Each additional selected candidate minimizes 
the sum of the squares of its similarities to all prior selected candidates. 


SIMULATION RESULTS 

For the initial tests of the spectral correction algorithm, the model 
shortwave spectral radiances are represented as the corresponding spectral 
radiant exitances reduced by the Lambertian factor tt. Instrument errors 
other than spectral dependence are not considered. 
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Simulation of Clear Sky Mid-Latitude Ocean Radiance Estimation 

The Earth radiation fields consist of six shortwave fields and nine 
longwave fields which are consistent with clear mid-latitude ocean scenes. 

Each shortwave-longwave field pair is assigned the same probability (1/54). 

For the shortwave fields, the model atmospheres used in the radiation transfer 
calculations are Wark's (ref. 3) models 4, 29, 34, 42, 44, and 45. For the 
longwave fields, we use models 2, 4, 29, 32, 34, 42, 44, 45, and 101. The 
shortwave fields assume a solar zenith angle of 45°. 


Table l.a shows the results of 25 pseudo-random samples of the clear mid- 
latitude scene class. The radiance estimate errors are all very small, less 
than the expected instrument calibration errors. Both the sample and popula- 
tion error statistics are summarized in Table 1 .b. The influence matrix (the 
B matrix of equation (16)), also displayed in Table l.b, indicates that the 
longwave channel is given little influence on the radiance estimates. This is 
understandable in view of the longwave channel's relatively large variation in 
response to unit longwave radiance, with a population standard deviation of 
0.02247, which is 3.9 percent of the mean longwave response, an order of 
magnitude more than that of the other channels for either shortwave or long- 
wave fields. 

The nighttime case is simulated by nulling the shortwave fields. The 
estimation errors (Tables 2. a and 2.b) are typically smaller than the daylight 
errors. Again the longwave channel has small influence. 


Simulation of Ambiguous Clear/Hazy Mid— Latitude Ocean Radiance Estimation 

The Earth radiation fields consist of the ones used for the clear sky 
case plus three shortwave fields consistent with haze over mid-latitude ocean 
and two longwave fields consistent with mid-latitude ocean covered by clouds 
at or lower than the 800-mbar level. For the additional longwave fields, the 
model atmospheres used in the radiation transfer calculations are Wark's 
(ref. 3) models 83 and 97, with cloud decks imposed at the 800- and 850-mbar 
levels, respectively. For the additional shortwave fields, we use models 2, 
32, and 101. The shortwave fields assume a solar zenith angle of 45°. Each 
shortwave-longwave field pair is assigned the same probability. 

The daytime case is summarized n Tables 3. a and 3.b. Again the estima- 
tion errors are smaller than the expected instrument calibration errors. The 
elements of the influence matrix are moderately different from those of the 
previous daylight case. 

The nighttime case (Tables 4.a and 4.b) longwave estimation errors are 
comparable to the daylight longwave errors and to the nighttime clear ocean 
longwave errors. The influence matrix is nearly identical to the nighttime 
clear ocean matrix. 
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CONCLUDING REMARKS 


Simulations of the spectral correction algorithm have shown it to be 
stable and to yield negligible errors when the scene spectral fields are well 
characterized* The limitations of the method derive from uncertainties in the 
spectral response of the channels and any lack of realism in the model spec- 
tral fields and their associated probabilities. The spectral response is of 
particular concern because the instruments must be characterized spectrally 
before flight; the in-flight calibration is not sufficient to accomplish this* 

< simulation studies presented are highly idealized^ but they serve to 
confirm the soundness of the algorithm. 
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Table l.a 


]PR*Q 

ASSUMED FRACTIONAL TRUNCATION ERROR IN YY ELEMENTS^ IE-08 
THE FIELD PROBABILITY MATRIX FILE NAME IS I, 0CM.Z45.PR0B2<1 , ‘ 

THE SCANNER RESPONSE MATRIX FILE NAME IS ,, MID-LAT.OCEAN.36FIELDS' 1 
THERE ARE 25 MEASUREMENT TRIPLETS. 


ACTUAL SW 

EST SW 

SW ERROR 

ACTUAL 

LW EST LW 

LW ERROR 

29.633 

29.686 

-.052 

77.040 

77.068 

-.028 

29.777 

29.815 

-.037 

60.437 

60.467 

-.030 

29.641 

29.688 

-.046 . 

80.417 

80.457 

-.039 

29.633 

29.686 

-.052 

77.040 

77.068 

-.028 

29.777 

29.816 

-.038 

77.040 

77.070 

-.030 

29.330 

29.358 

-.027 

60.532 

. 60.556 

-.023 

29.777 

29.812 

-.034 

80.417 

80.459 

-.041 

29.562 

29.509 

.053 

85.020 

85.049 

-.029 

29.777 

29.815 

-.037 

60.532 

60.562 

-.029 

29.562 

29.510 

.051 

97.985 

98.051 

-.065 

29.633 

29.684 

. -.051 

60.437 

60.465 

-.027 

29.641 

29.692 

-.050 

60.532 

60.560 

-.027 

29.562 

29.504 

.057 

60.437 

60.468 

-.030 

29.777 

29.809 

-.031 

54.147 

54.191 

-.043 

29.633 

29.689 

-.055 

85.020 

85.046 

-.026 

29.330 

29.358 

-.028 

71.543 

71.564 

-.021 

29.777 

29.809 

-.031 

54.147 

54.191 

-.043 

29.330 

29.364 

-.033 

97.985 

98.044 

-.058 

29.641 

29.695 

-.053 

87.532 

87.571 

-.039 

29.777 

29.821 

-.043 

97.985 

98.050 

-.064 

29.330 

29.360 

-.030 

87.532 

87.567 

-.035 

29.330 

29.358 

-.027 

60.532 

60.556 

-.023 

29.562 

29.505 

.056 

77.040 

77.071 

-.031 

29.777 

29.819 

-.042 

85.020 

85.049 

-.028 

29.633 

29.687 

-.054 

87.532 

87.571 

-.039 


/ 



Table l.b 


THE SW RADIANCES ARE BIASED BY .025 OR .086% 

(PERCENTAGES ARE OF MEAN SW OR LW RADIANCES.) 

THE LW RADIANCES ARE BIASED BY .035 OR .047 % 

THE SW STANDARD DEVIATION IS .044 OR .150 % 

THE LW STANDARD DEVIATION IS .037 OR .049 % 

THE SW STANDARD DEVIATION EXCLUDING BIAS IS .036 

THE LW STANDARD DEVIATION EXCLUDING BIAS IS .012 


MATRIX OF CHANNEL INFLUENCES ON ESTIMATES 


TOTAL CHANNEL SW CHANNEL 


LW CHANNEL 


SW EST 1.05237961E-03 1.6670959 

LW EST 1.19022346 -1.35070884 


- 1.204 9973E-03 
-.127950311 


SW POPULATION BIAS* .01389 
LW POPULATION BIAS* .03614 
SW POPULATION STD DEV= .04428 
LW POPULATION STD DEV= .03800 


OR .046 % 

OR .048 % 

OR .149 % 
OR .050 % 



Table 2. a 


ASSUMED FRACTIONAL TRUNCATION ERROR IN YY ELEMENTS= IE-08 
3 BY 3 YY MATRIX IS SINGULAR. 

RADIANCE ESTIMATES ARE DERIVED FROM LW CHANNEL AND TOTAL CHANNEL. 
THE FIELD PROBABILITY MATRIX FILE NAME IS ,, OCM.Z45.PROB24" 

THE SCANNER RESPONSE MATRIX FILE NAME IS “MID-LAT. OCEAN. 36FIELDS" 
THERE ARE 25 MEASUREMENT TRIPLETS. 


ACTUAL SV 

EST SV 

SV ERROR 

ACTUAL LV 

EST LV 

LV ERROR 

0 

0 

0 

54.147 

54.150 

-2.716E-03 

0 

0 

0 

77.040 

77.029 

.010 

0 

0 

0 

80.417 

80.385 

.032 

0 

0 

0 

80.417 

80.385 

.032 

0 

0 

0 

97.985 

98.022 

-.037 

0 

0 

0 

71.543 

71.539 

4.287E-03 

0 

0 

0 

54.147 

54.150 

-2.716E-03 

0 

0 

0 

87.532 

87.530 

1.544E-03 

0 

0 

0 

60.532 

60.538 

-5.268E-03 

0 

0 

0 

60.437 

60.442 

-5.050E-03 

0 

0 

0 

97.985 

98.022 

-.037 

0 

0 

0 

85.020 

85.015 

5.385E-03 

0 

0 

0 

77.040 

77.029 

.010 

0 

0 

0 

80.417 

80.385 

.032 

0 

0 

0 

97.985 

98.022 

-.037 

0 

0 

0 

60.532 

60.538 

-5.268E-03 

0 

0 

0 

80.417 

80.385 

.032 

0 

0 

0 

97.985 

98.022 

-.037 

0 

0 

0 

60.437 

60.442 

-5.050E-03 

0 

0 

0 

80.417 

80.385 

.032 

0 

0 

0 

80.417 

80.385 

.032 

0 

0 

0 

87.532 

87.530 

1.544E-03 

0 

0 

0 

77.040 

77.029 

.010 

0 

0 

0 

71.543 

71.539 

4.287E-03 

0 

0 

0 

60.437 

60.442 

r 5.050E-03 



• Table 2. b 


(PERCENTAGES ARE OF MEAN £W OR LW RADIANCES.) 

THE LW RADIANCES ARE BIASED BY — 2.629E-03 OR -3.425E-03 % 

THE LW STANDARD DEVIATION IS .022 OR .029 % 

THE LW STANDARD DEVIATION EXCLUDING BIAS IS .022 

MATRIX OF CHANNEL INFLUENCES ON ESTIMATES 

TOTAL CHANNEL SW CHANNEL LW CHANNEL 

SW EST 0 0 0 

LW EST 1.22436789 0 —.174493093 

SW POPULATION BIAS= 0 

LW POPULATION BIAS= -5.02077E-04 OR -6.69776E-04 % 

SW POPULATION STD DEV= 0 

LW POPULATION STD DEV= .01719 OR .02294 % 


r 



Table 3. a 


ASSUMED FRACTIONAL TRUNCATION ERROR IN YY ELEMENTS= IE-08 
THE FIELD PROBABILITY MATRIX FILE NAME IS , 'OCHM.Z45.PROB27‘' 

THE SCANNER RESPONSE MATRIX FILE NAME IS ' ‘MID -LAT. OCEAN. 3 6FI ELDS" 
THERE ARE 25 MEASUREMENT TRIPLETS. 


ACTUAL SW 

EST SW 

SW ERROR 

ACTUAL 

LW EST LW 

LW ERROR 

29.562 

29.492 

.070 

87.532 

87.535 

-3 . 829E-03 

29.330 

29.360 

- .029 

97.985 

98.021 

-.036 

29.330 

29.322 

8.158E-03 

80.417 

80.389 

.028 

29.562 

29.492 

.070 

78.784 

78.779 

5.619E-03 

37.171 

37.117 

.054 

77.040 

77.040 

-2.762E-04 

37.171 

37.108 

.063 

54.147 

54.163 

-.015 

29.777 

29.814 

-.036 

97.985 

98.022 

-.037 

37.691 

37.768 

-.077 

87.532 

87.536 

-4.525E-03 

29.777 

29.800 

-.022 

78.784 

78.774 

.010 

37.691 

37.767 

-.075 

77.040 

77.035 

5.330E-03 

37.275 

37.262 

.012 

60.437 

60.449 

-.011 

37.691 

37.744 

-.052 

80.417 

80.395 

.022 

29.641 

29.667 

-.026 

54.147 

54.151 

-3.738E-03 

29.633 

29.671 

-.037 

78.784 

78.773 

.010 

29.562 

29.492 

.070 

87.532 

87.535 

-3.829E-03 

29.641 

29.684 

-.042 

85.020 

85.013 

7.351E-03 

29.562 

29.506 

.056 

97.985 

98.026 

-.041 

29.777 

29.800 

-.023 

87.532 

87.531 

6.019E-04 

37.275 

37.262 

.012 

60.437 

60.449 

-.011 

29.330 

29.350 

-.019 

65.705 

65.702 

3. 138E-03 

29.562 

29.495 

.066 

65.705 

65.707 

-2.281E-03 

37.691 

37.768 

-.076 

78.784 

78.779 

4 923E-03 

29.777 

29.807 

-.029 

85.020 

85.014 

6.236E-03 

29.330 

29.322 

8. 158E-03 

80.417 

80.389 

.028 

29.330 

29.346 

-.016 

87.532 

87.530 

1.590E-03 


f 



Table 3.b 


THE SV RADIANCES ARE BIASED BY 2.937E-03 OR 9.156E-03V. 

(PERCENTAGES ARE OF MEAN SW OR LW RADIANCES.) 

THE LW RADIANCES ARE BIASED BY 1.490E-03 OR 1.888E-03 % 

THE SW STANDARD DEVIATION IS .048 OR .151 % 

THE LW STANDARD DEVIATION IS .017 OR .022 % 

THE SW STANDARD DEVIATION EXCLUDING BIAS IS .048 

THE LW STANDARD DEVIATION EXCLUDING BIAS IS .017 


MATRIX OF CHANNEL INFLUENCES ON ESTIMATES 


TOTAL CHANNEL SW CHANNEL 


LW CHANNEL 


w SW EST .0224043801 

LW EST 1.22373034 


1.63077869 

-1.40781623 


-.0300852433 

-.173524708 


SW POPULATION BIAS= 7.61532E-04 
LW POPULATION BIAS= 2.04801E-03 
SW POPULATION STD DEV= .04731 
LW POPULATION STD DEV= .01564 


OR 2.365E-03 % 
OR 2.750E-03 % 
OR .146 % 
OR .021 'A 


t < 


f ’ 



Table 4. a 


ASSUMED FRACTIONAL TRUNCATION ERROR IN YY ELEMENT£= IE-08 
3 BY 3 YY MATRIX IS SINGULAR. 

RADIANCE ESTIMATES ARE DERIVED FROM LV CHANNEL AND TOTAL CHANNEL. 
THE FIELD PROBABILITY MATRIX FILE NAME IS M OCHM. 245^0327“ 

THE SCANNER RESPONSE MATRIX FILE NAME IS “MID-LAT. OCEAN. 36FIELDS" 
THERE ARE 25 MEASUREMENT TRIPLETS. 


ACTUAL 5V 

EST SV 

SV ERROR 

ACTUAL LV 

EST LV 

LV ERROR 

0 

0 

0 

65.705 

65.704 

7. 131E-04 

0 

0 

0 

65.705 

65.704 

7.131E-04 

0 

0 

0 

71.543 

71.540 

2.979E-03 

0 

0 

0 

60.532 

60.539 

-6.413E-03 

0 

0 

0 

87.532 

87.532 

7.748E-07 

0 

0 

0 

54.147 

54.151 

-3.754E-03 

0 

0 

0 

54.147 

54.151 

-3.754E-03 

0 

0 

0 

85.020 

85.016 

3.873E-03 

0 

0 

0 

78.784 

78.774 

.010 

0 

0 

0 

78.784 

78.774 

.010 

0 

0 

0 

54.147 

54.151 

-3.754E-03 

0 

0 

0 

65.705 

65.704 

7. 131E-04 

0 

0 

0 

80.417 

80.386 

.031 

0 

0 

0 

97.985 

98.024 

-.038 

0 

0 

0 

71.543 

71.540 

2. 979E-03 

0 

0 

0 

97.985 

98.024 

-.038 

0 

0 

0 

54.147 

54.151 

-3. 754E-03 

0 

0 

0 

71.543 

71.540 

2.979E-03 

0 

0 

0 

60.437 

60.443 

-6. 193E-03 

0 

0 

0 

60.532 

60.539 

-6.413E-03 

0 

0 

0 

87.532 

87.532 

7.748E-07 

0 

0 

0 

65.705 

65.704 

7.131E-04 

0 

0 

0 

77.040 

77.031 

9.376E-03 

0 

0 

0 

71.543 

71.540 

2.979E-03 

0 

0 

0 

85.020 

85.016 

f .3.873E-03 



Table 4.b 


(PERCENTAGES ARE OF MEAN SW OR 
THE LW RADIANCES ARE BIASED BY 1.095E- 
THE LW STANDARD DEVIATION IS .013 


LW RADIANCES.) 

03 OR 1.519E-03 S 

OR .oie % 


THE LW STANDARD DEVIATION EXCLUDING BIAS IS .013 


MATRIX OF CHANNEL INFLUENCES ON ESTIMATES 


TOTAL CHANNEL 

SW EST 0 

LW EST 1.22441104 


SW CHANNEL 

0 

0 


LW CHANNEL 

0 

-.174529083 


ro 


SW POPULATION BIA£= 0 

LW POPULATION BIAS= -3.23628E-04 OR -4.34586E-04 % 

SW POPULATION STD DEV= 0 

LW POPULATION STD DEV= .01593 OR .02139 % 
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